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Ultrahigh Energy Cosmic Rays and Neutrinos 
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We discuss the relation between the highest energy cosmic rays (UHECR) and UHE neutrinos. The neutrinos 
produced in the sources of optically thin astrophysical sources have been linked to the UHECR emissivity of the 
Universe. The fluxes of cosmogenic neutrinos, generated in propagation by UHECR, also reflect the acceleration 
of these particles, the maximum acceleration energy, and the cosmological evolution of their sources. 



1. Introduction 

It is now a time when both subjects in the ti- 
tle are of high scientific interest in the field of 
particle astrophysics. The interest is related to 
the emergence of new set of experiments that are 
obviously of much better quality than what we 
know from the past. 

The Southern Auger Observatory [Tl in Ar- 
gentina is almost fully completed and is operat- 
ing during construction. The Telescope Array [2] 
(TA) is under construction in Utah. Both these 
giant air shower arrays are designed and con- 
structed to operate in hybrid mode, i.e. they em- 
ploy both methods of air shower detection: sur- 
face array plus optical detectors that follow the 
fluorescent light generated by the giant air shower 
in the atmosphere. The optical detectors operate 
about 10% of the time but the energy assignment 
by integration over the shower longitudinal profile 
seems to be less model dependent than the classi- 
cal method of relating the primary energy to the 
shower particle density on the ground. The latter 
depends much more on the hadronic interaction 
model used in Monte Carlo studies of the shower 
development. The fraction of showers detected in 
hybrid mode by both detectors can be used for a 
better calibration of the surface array. The sizes 
of the surface arrays, 3,000 km^ for Auger and 
1,000 km^ for TA are very impressive, more than 
one order of magnitude higher than the previous 
largest array - Agasa. 

There are also big news on the neutrino front. 



The IceCube neutrino observatory [3] on the 
South Pole is in the middle of successful construc- 
tion. Twenty two of the eighty underice strings 
are deployed and the deployment is now faster 
than initially expected. The Antares [4] detector 
in the Mediterranean is also successfully close to 
completion and the kmSnet [5] is working on the 
design of a cubic kilometer under water neutrino 
telescope in the Northern hemisphere. With these 
cubic kilometer neutrino detectors we approach 
for first time the dimensions needed for detection 
of astrophysical neutrinos - a gigaton of target 
matter or 6xl0'^* target nucleons. 

Experiments of this size and quality can cer- 
tainly make the objects of their studies fashion- 
able. The first high statistic results of the Auger 
observatory were posted on ArXive.org and re- 
ported at the 20th International Cosmic Ray Con- 
ference in Merida, Mexico in front of a huge au- 
dience. Auger has set [?] a strict limit of 2% 
on the gamma-ray contribution to the UHECR 
flux above 10^^ eV and in the following we shall 
assume that the highest energy cosmic rays origi- 
nate in acceleration in powerful astrophysical ob- 
jects. 

IceCube has not yet detected any astrophysical 
neutrino signals but the limits already obtained 
have let to revisions of many models. 



2. Relation and differences between UHE 
cosmic rays and neutrinos 

The connection between ultrahigh energy cos- 
mic rays (UHECR) and neutrinos was first em- 
phasized by Waxman&Bahcall [7]- This paper 
used an estimate of the cosmic ray emissivity of 
the Universe :8, above 10^^ eV to calculate the 
maximum neutrino fluxes generated at the accel- 
eration sites of UHECR. The calculation has to 
assume an energy spectrum for these particles. 
It was taken to be a flat E^'^ acceleration spec- 
trum. Accounting for the cosmological evolution 
of the UHECR sources (and assuming that the 
same acceleration spectrum is followed to much 
lower energy) the maximum isotropic neutrino 
flux was calculated to be E'^dF/dE = 5x 10~^ per 
cm^.ster.s and was called the upper bound of the 
isotropic neutrino fluxes. This bound applies to 
optically thin sources where the accelerated nuclei 
can leave the astrophysical source and does not 
limit the neutrino production in optically thick 
ones. 

A better calculation emphasizing the sim- 
plifications in this approach was published by 
Mannheim, Protheroe & Rachen [9 . Their up- 
per hound touches the W&B one at E^ = 10^* eV 
and is higher than it at all other neutrino ener- 
gies. A comparison between the two calculations 
is shown in Fig. [TJ Neither of these calculation 




Figure 1 . The upper bounds of astrophysical neu- 
trinos calculated by Waxman&Bahcall (straight 
line) and Mannheim, Protheroe & Rachen. 



include the neutrinos that are generated by the 
UHECR on propagation to us - the cosmogenic 
neutrinos. 

Ultrahigh energy cosmic rays are not well un- 
derstood for two reasons: it is indeed very diffi- 
cult to appreciate how charged particles can be 
accelerated to energies as high as 10^° eV in any 
astrophysical object, and because of their high en- 
ergy loss in interactions on the microwave back- 
ground (MBR) and other universal photon fields, 
such as the infrared/optical background (IRB). 
Although we are discussing interactions of par- 
ticles exceeding 10^^ eV in the Lab the physics 
involved is very well known. Because of the low 
energy of the photon background fields the cen- 
ter of mass energy is in the 1-10 GeV range and 
the cross sections are well known from accelerator 
experiments. There are two important processes: 
photoproduction interactions and e+e^ pair pro- 
duction. The photoproduction process that gen- 
erates the GZK [in] effect has in the MBR a 
threshold of about 3x10^^ eV for protons. There 
are lower energy proton interactions in the IRB 
but they do not seem to be important for the 
proton spectrum evolution in propagation. The 
minimum interaction length is at proton energy 
of 5x10^° eV and is below 4 Mpc. At higher 
energy the energy loss length is about 14 Mpc, 
which for Hubble constant hg = 0.75 corresponds 
to redshift z=0.0035. A large fraction of the pro- 
ton energy loss in photoproduction interactions 
goes into neutrinos. Every time a charged pion 
is produced three neutrinos (muon neutrino and 
antineutrino and an electron neutrino) arc pro- 
duced. 

At lower energy (smaller CM energy needed) 
the main energy loss process is the electron- 
positron pair production. This process has a 
threshold of 2x10^^ eV. The cross section in- 
creases with proton energy, but the energy loss 
per interaction decreases. The combination 
leads to a minimum energy loss distance of 1.2 
Gpc (z = 0.3) at about 2x10^9 eV. Berezin- 
sky&Grigorieva [11] first discussed the impor- 
tance of this process for the UHECR spectrum 
after propagation. 

At energies lower than 2x10^* eV the main en- 
ergy loss process is the adiabatic energy loss from 



the expansion of the Universe. The energy loss 
length for /iq = 0.75 is 4 Gpc. 

Both the UHECR spectra from isotropically 
distributed cosmic ray sources and the neutrinos 
generated by them depend on the cosmic ray ac- 
celeration spectrum and on the maximum accel- 
eration energy. There is, however, a big difference 
related to the very different energy loss of cosmic 
rays and neutrinos. Neutrinos only suffer adia- 
batic energy loss and easily propagate to us from 
all redshifts while protons of arbitrary energy can 
propagate from large redshifts. Figure [2] shows 
the proton spectrum after propagation at differ- 
ent redshifts, from 0.025 to 0.4 The maximum 
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Figure 2. Evolution of the cosmic ray spectrum 
in propagation on different distances. 

energy in this graph is 10^^ eV with exponential 
cutoff at 10^^'^ eV. Still, after propagation on z = 
0.4, no protons of energy above 10^^ remain in the 
cosmic ray flux. This fact shows that the influ- 
ence of the cosmological evolution of the UHECR 
sources is not very significant for the highest en- 
ergy cosmic rays. 

The cosmogenic neutrino fluxes, on the other 
hand, are very sensitive to the cosmological evo- 
lution of the sources. Figure [3] shows a similar 
graph for cosmogenic neutrinos. The contribu- 
tion of different redshifts is shown on a logarith- 
mic scale of z for the same cosmological evolu- 
tion model as in Fig. [2]- (1 + z)'^ to z = 1.9 and 
then constant to 2.8 with exponential decline at 



higher redshift. Because of the small energy loss 
of neutrinos the contribution continues growing 
after redshift of 2.5. Without cosmological evolu- 
tion of the sources the highest contribution would 
have come from the contemporary Universe. If 




Figure 3. Contribution of different redshifts to 
the flux of cosmogenic neutrinos. One can see 
the adiabatic energy loss of neutrinos as the peak 
of the distributions change with redshift. 



we succeed in the detection of cosmogenic neutri- 
nos we can understand the cosmological evolution 
of the extragalactic cosmic ray sources and cre- 
ate a better model of the highest energy cosmic 
rays [12]. 

3. Ultrahigh Energy Cosmic Rays 

In this section we will present the newest set of 
data, compare it to older ones and to some of the 
available models. Figure [5] shows the data from 
all experiments. The new Auger data come from 
the surface array normalized to the fluorescent 
detector energy assignment in hybrid events |13j . 
At lower energy only showers detected in hybrid 
mode are included [14] . Another energy spec- 
trum, consistent with the shown ones, was derived 
from inclined showers. At 10^^ eV the difference 
in the energy assignment between the highest and 
lowest flux (Agasa [15^ and Auger) is about 40%. 
Auger supports the measurement of HiRes |16ll7j 
that shows a strong decrease of the cosmic ray 



flux above 5-6x10^® eV. The total number of 
events higher than 10^° eV is two. These two 
data sets change our expectations for such events 
- it appears that expect to see 0.5 event in 1,000 
km^.ster.yr of exposure, about 10 times less than 
the Agasa estimate. The Auger energy assign- 
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Figure 4. Spectrum of the highest energy cosmic 
rays detected by different experiments. 

ment seems to be somewhat lower than that of 
HiRes, maybe by about 20% averaged over the 
whole energy range shown. There are also some 
minor differences in the exact shape of the spec- 
trum: the dip at lO^*'^ is better pronounced and 
the recovery at higher energy is faster. The GZK 
effect seems to take place at slightly lower en- 
ergy and the flux decrease is not as steep as in 
HiRes data set. As far as the highest energies are 
concerned the shower statistics is too small to be 
analyzed. 

The question now is: which of the available 
models fit the experimental spectrum the best. 
There are in principle three available models. The 
one of Berezinsky et al. [18I19J uses steep accel- 
eration spectrum (E~^-^) down to about 10"'^^ eV. 
This model emphasizes the dip in the spectrum 
that is due to the e+e^ pair production loss in 
propagation and its conversion to purely adia- 
batic energy loss. This model is unique because it 
fits the cosmic ray flux above 10^^ eV with extra- 



galactic cosmic ray protons without any galactic 
component. The model fits the HiRes spectrum 
better than it fits the Auger one. 

The second, also proton model, was suggested 
by Waxman&Bahcall [20 and supported also by 
other authors. It uses a fiat E~^ acceleration 
spectrum, strong cosmological evolution of the 
cosmic ray sources and predicts a dip where the 
flux of extragalactic cosmic rays intersects the 
galactic cosmic ray component. This model re- 
quires that our Galaxy accelerates some cosmic 
rays to energies above 10^^ eV. 

A third model uses cosmic rays at their sources 
with a composition similar to the GeV galac- 
tic cosmic rays |21I22| . Since there is a signifi- 
cant fraction of heavy cosmic ray nuclei the en- 
ergy loss in propagation is different in this model. 
Heavy nuclei lose energy in photodisintegration 
in the photon fields. The energy threshold co- 
incides with the giant dipole resonance, and is 
thus much lower than photoproduction. The nu- 
clei lose nucleons in the process - the energy per 
nucleon is stable but the total energy per nucleus 
decreases. The total energy loss length is not dis- 
similar to that for proton photoproduction. As 
well as the other two models it fits the measured 
spectra quite well with an intermediate accelera- 
tion spectrum of E~^-^~^'^. 

The three models predict quite different nuclear 
compositions for the UHECR. In the Berezinsky 
et al model the transition from galactic to extra- 
galactic cosmic rays, and from very heavy to very 
light composition, is at the approach of 10^* eV. 
In the flat spectrum model this transition is at 
significantly higher energy and UHECR should 
contain some iron nuclei even above 10^^ eV. The 
mixed composition model is somewhat intermedi- 
ate. In the whole energy range there is a signifi- 
cant fraction of protons released in the photodis- 
integration process complementing the primary 
and secondary nuclei. Only in the highest energy 
range the composition becomes light. 

The experimental data shown in Fig. [5] do 
not seem to support any of the models. The 
circles and the pentagons show the HiRes and 
HiRes/MIA data [23] on the shower depth of 
maximum converted to {In A) using the Siby 112.1 
hadronic interaction model. The diamonds show 
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Figure 5. Cosmic ray nuclear composition at the 
highest energies. The measured quantity is the 
depth of maximum which is converted by me to 

{In A). 



the the Auger data [25 treated the same way. 
The composition measurement of Auger shows 
generally a heavier composition, consistent with 
Helium. Excluding the highest energy points 
the differences between Auger and HiRes are not 
huge. The beautiful picture presented in Ref. [23] 
does not exist any more, and the composition pic- 
ture seems to be much more complicated. 

The composition studies done by Auger (and 
probably HiRes) are hurt by the lack of under- 
standing of the air shower development revealed 
by the hybrid detection method. If Auger did 
not have a fluorescent detector and used the sur- 
face shower densities as Agasa did, they would 
be two changes in the current results: the energy 
estimate would be higher by about 25% and the 
cosmic ray composition would appear much heav- 
ier [25j . The currently used hadronic interaction 
models are not able to explain this effect. 

One important result from the studies of the 
composition of UHECR is the strict limit of the 
fraction of 7-rays which is set to 2% for all par- 
ticles above 10^^ eV f5]. This limit shows that 
at least 98% of the UHECR particles are nuclei 
accelerated in astrophysical objects. 



4. Cosmogenic neutrinos 

Cosmogenic neutrinos are the neutrinos gener- 
ated in photoproduction interactions of the prop- 
agating cosmic rays in the photon fields of the 
Universe. In the MBR the current threshold en- 
ergy for proton photoproduction is about 3x10^^ 
eV. Protons lose less than 20% of their energy 
in the threshold energy range. We can roughly 
estimate the average cosmogenic neutrino energy 
as {El,) = Ep X K„iei/4: where K^nei is the av- 
erage energy loss of the protons per interaction 
and each neutrino takes 1/4 of the pion energy, 
i.e. {Ep) — 10^'* eV. Cosmogenic neutrinos were 
first proposed by Berezinsky & Zatsepin [26j and 
have been since calculated many times, most re- 
cently in Ref. [27] . Neutrinos generated at higher 
redshift have adiabatic losses. In addition the cos- 
mological evolution of MBR makes possible the 
interactions of lower energy protons so the aver- 
age energy of the cosmogenic neutrinos after in- 
tegration in redshift is of the same order. 

The spectrum of cosmogenic neutrinos de- 
pends on the UHECR acceleration spectrum, the 
UHECR source distribution and very strongly 
on the cosmological evolution of the UHECR 
sources \TT . Flat acceleration (7=1) models gen- 
erate high fiux because they contain higher num- 
ber of interacting protons of energy above 3x 10^^ 
eV and because they need strong cosmological 
evolution of the cosmic ray sources in order to 
fit the observed cosmic ray spectrum. The model 
of Berezinsky et al, on the other hand, has much 
smaller number of interacting protons for the 
same source emissivity and does not need cosmo- 
logical evolution of the sources. Figure [6] shows 
the fluxes of cosmogenic neutrinos generated by 
these models. Note that the De spectrum peaks 
at about 10^^'^ eV while all other flavors peak 
at 10^^ eV. The reason is that D^. are generated 
in neutron decay rather than in photoproduction 
interaction. They take a very small fraction of 
the neutron energy. Since the acceleration spec- 
trum is protons only, the neutrons are secondaries 
coming from charge exchange interactions. Only 
at very high energy some secondary neutrons in- 
teract to produce high energy Vg with the same 
distribution as the other flavors. 
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Figure 6. Cosmogenic neutrinos generated by 
protons with a flat acceleration spectrum (7=1.0) 
with cosmological evolution and a steep one 
(7=1.7) without evolution. Electron neutrinos 
are shown with a solid line, electron antineutrinos 
- with dots, muon neutrinos - with dashes and 
muon antineutrinos with dash-dotted line. The 
symbols show the sum of all neutrino flavors. 



Figure 7. Cosmogenic neutrinos generated by 
protons with a flat acceleration spectrum (7=1.0) 
with cosmological evolution and a steep one 
(7=1.7) without evolution in MBR & IRB. The 



pentagons show the spectrum of v^ 



generated 



by the mixed composition model and the x's show 
the number of v^. generated by the same model. 



It is of some importance to note that MBR is 
not the only target for neutrino production. The 
second most important one is the isotropic in- 
frared and optical background (IRB). Its number 
density is, or course, much lower, but lower en- 
ergy protons can interact in it and even in the 
case of flat acceleration spectra the number of in- 
teracting protons to a large extent compensates 
for the lower photon target density. 

Fig. [7] shows the spectra of cosmogenic neu- 
trinos generated by a flat (7=1.0) and steep 
(7=1.7) UHECR acceleration spectra in the MBR 
and IRB. Since the steep injection spectrum has 
higher number of lower energy protons it provides 
more interactions in the IRB and decreases the 
difference between the two models. The flat ac- 
celeration spectrum model contains more high en- 
ergy cosmogenic neutrinos 

The mixed composition model generates 
mostly Ve due to the decaying neutrons released 
by the photodisintegrating neutrons [28]. Since 
there are many more such neutrons than pho- 
toproduction interactions electron antineutrinos 
dominate the cosmogenic neutrino flux. Electron 



neutrinos would have the same spectrum as v^ 
and t?^, the peak of which is similar to the neutri- 
nos from the flat spectrum model. Note that such 
much lower energy neutrinos have a significantly 
lower cross section and do not contribute much to 
the total event rate. Their peak energy is also be- 
low the the Glashow resonance energy of 6x10^ 
GeV. Still, in case we are lucky enough to de- 
tect cosmogenic neutrinos with the neutrino tele- 
scopes under construction and design they could 
help a lot in limiting the models for the origin 
of the ultrahigh energy cosmic rays. We have 
to remember, however, the expected event rate 
of cosmogenic neutrinos is small, less than 1 per 
km'^.yr and such detection requires new detector 
technologies (radio detection?) that can cover 
hundred km^. 

5. Summary 

The new UHECR data sets agree that UHECR 
spectrum experiences a steepening energy spec- 
trum at about 5-6x10^^ eV which looks consis- 
tent with the GZK effect which is due to the 
energy loss of these particles in photoproduction 



or photodisintegration interactions. It is not yet 
obvious what are the parameters of the acceler- 
ation process and what is the evolution of the 
cosmic ray sources. All these questions will prob- 
ably have to wait until the total world statistics 
is increased by a large factor. 

The composition of all particles above 10^^ eV 
is not yet established. Different experimental 
data are not far away from each other, but do 
not currently contribute to the understanding of 
the UHECR origin. The advance in that respect 
is the strict limit of 2% that the Auger collabora- 
tion has set on the fraction of 7-rays in the cosmic 
ray flux above 10^^ eV. 

Different UHECR models predict various fluxes 
of cosmogenic neutrinos that are generated by the 
cosmic rays at propagation from their sources to 
us. Detection of such neutrinos and a compar- 
ison of their fluxes to the direct observations of 
UHECR would contribute significantly to under- 
standing of the UHECR origin. For this to hap- 
pen, however, we will have to rely on new, bigger 
neutrino observatories. 
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